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ABSTRACT 



Context. There is a long-standing controversy concerning the reason why the Mn i 5394.7 A line in the solar irradiance spectrum 
brightens more at larger activity than most other photospheric lines. The claim that this activity sensitivity is caused by spectral inter- 
locking to chromospheric emission in Mg n h & k is disputed. 
Aims. To settle and close this debate. 

Methods. Classical one-dimensional modeling is used for demonstration; modern three-dimensional MHD simulation for verification 
and analysis. 

Results. The Mm 5394.7 A line thanks its unusual sensitivity to solar activity to its hyperfine structure. This overrides the thermal 
and granular Doppler smearing through which the other, narrower, photospheric lines lose such sensitivity. We take the nearby Fe i 
5395.2 A line as example of the latter and analyze the formation of both lines in detail to demonstrate and explain granular Doppler 
brightening. We show that this affects all narrow lines. Neither the chromosphere nor Mg n h & k play a role, nor is it correct to 
describe the activity sensitivity of Mn i 5394.7 A through plage models with outward increasing temperature contrast. 
Conclusions. The Mm 5394.7 A line represents a proxy diagnostic of strong-field magnetic concentrations in the deep solar pho- 
tosphere comparable to the G band and the blue wing of Ha, but not a better one than these. The Mn i lines are more promising as 
diagnostic of weak fields in high-resolution Stokes polarimetry. 

Key words. Sun: photosphere - Sun: chromosphere - Sun: granulation - Sun: magnetic fields - Sun: faculae, plage 



1. Introduction 

In this paper we analyze the formation of the solar Mm 
5394.7 A line in the context of its global sensitivity to so- 
lar activity, a subject which has received considerable attention 
based on the extensive observations of t his line by Livingston 
and coworkers and Vince and coworkers (iLivingston & Wallace] 
| 1987l:IVince&ErkaDidll998l: Danilovic & Vince 2004, 2005; 
Malanushenko et al. 2004; Danilovic et al. Vince et al. 

2005a, 2005b: lLivingston et al.ll2007h . 

In hindsight, the starting publications were by iThackerayl 

(Il937l) who pointed out that the violet wing of Mg n k (line cen- 
ter k 3 at 2795.53 A) overlaps with Mm 2794.82 A and so may 
produce optical pumping of that and other Mn i lines in stellar 
spectra, and by lAbtl (Il952l) who pointed out that the unusual 
widths of the Mn i lines in the solar spectrum are due to hyperfine 
structure. These two points are key ones in almost all later work 
on solar Mm lines, effectively dividing this into two categories. 
The first addresses the unusual sensitivity of the Mm 5394.7 A 
line to global solar activity, mostly debating the forceful claim by 
iDovleetalJ (l2001h that this is explained by Thackeray's Mg n k 
coincidence operating in the solar chromosphere. This activity 
response is also our subject here, but we establish and explain 
that neither Mgn k nor the chromosphere has anything to do 
with it. 
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We first summarize the observations and discuss these con- 
flicting views in the next paragraphs of this introduction, and we 
then put the issue to rest by demonstrating that Abt's hyperfine 
structure is the key agent through reducing spectral-line sensi- 
tivity to thermal and convective Dopplershifts outside magnetic 
concentrations. 

The second category of solar Mm papers also addresses 
usage of Mm lines as diagnostic of solar magnetism but 
concentrates on quantitative measurement of weak internet- 
work fields exploiting the intricate line-center opacity varia- 
tion with wavelength imposed by hyperfine structure to disen- 
tangle weak- and strong-field signatures in Stokes polarimetry 
(Lopez Ariste et al. 2002, 2006a, 2006b; lAsensio Ramos et al] 
l2007USanchez Almeida et al.ll2008l) . This potentially more fruit- 
ful topic is not addressed here. 



1.1. Activity modulation 

Livingston's inclusion of the Mm 5394.7 A line into his long- 
term full-disk "sun-as-a-star" line profile monitoring from 1979 
onwards was suggested by Elste who had pointed out that their 
large hyperfine structure makes the Mm lines less sensitive to 
the questionable microturbulence parameter than other ground- 
state n eutral-metal lines that may serv e as temperature diag- 
nostic (lElste & Teskdll978l : lElstef [T987). Livingston found that 
this line is the only photospheric line in his full-disk monitor- 
ing that shows appreciable variation with the cycle, in good 
concert with the Can K full-disk intensity variation. Its equiv- 
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alent width in the irradianc e spectrum varies by up to 2% 
(Livi ngston & Wallace! fT987h . Figure 16 in the overview pa- 
per of iLivingston et al.l (120071) displays the variation in relative 
line depth (i.e., the minimum intensity of the line in the full- 
disk spectrum expressed as fraction of the continuum intensity 
and measured from the latter, plotted upside down and there- 
fore labeled "central intensity" in the caption). The same data 
are pl otted as relative line depth in Fig. 2 of iDanilovic et al.l 
(2007), overlaid by a theoretical modeling curve. The relative 
line-center intensity increases about 2% from cycle minimum 
to maximum, slightly more than the corresponding decrease in 
equivalent width. 

IVince et all (l2005al) used observations at the Crimea 
Observatory including Zeeman polarimetry to measure the 
changes in Mm 5394.7 A between plages with different appar- 
ent magnetic flux density. They found that the line weakens with 
inc reasing flux as concluded already by lElste & Teskd (Il978l) 
and lElstel (11987b . 

1.2. Chromospheric interpretation 

iDovle et al.l (l200lb gave their paper the title "Solar Mm 
5432/5395 A line formation explained' which we paraphrase in 
our title above. They used NLTE computations to claim that 
these Mm lines are sensitive to optical pumping through the 
overlap coincidence with Mg n k noted by Thackeray. The claim 
consisted of displaying Mm profiles for different solar atmo- 
sphere models without and with taking Mg n h & k into account. 
The different models specified ad-hoc variation in the onset 
heights of the chromospheric and transition-region temperature 
rises. Appreciable variation of the two Mn i lines was found and 
attributed to the spectral interlocking. 

However, closer inspection undermines this claim. Large 
changes, of order 30% in line depth, were shown to occur when 
Mg ii h & k and all other blanketing lines are not taken into ac- 
count, but this is not a relevant test since deletion of the ul- 
traviolet line haze invalidates the ionization equilibrium evalu- 
ation for any species with intermediate ionization energy. The 
same test would be as dramatic for any optical Fei line. The 
changes in the profile of Mm 5394.7 A between the cases of 
Mn i - Mg ii coupling and no coupling were negligible unless the 
model possessed a very deep-lying chromosphere and transition 
region producing unrealistic high peaks in Ca n H & K and Mg n 
h&k. Even then, the computed brightening of Mm 5394.7 A 
amounted to only a few percent, yet to be diluted through a fill- 
ing factor of order 0.01 to represent the contribution of active- 
sun plage in full-disk averaging. 

iDoyle et all d200l|) added no further analysis (such as spec- 
ification of NLTE departures, radiation fields, source functions) 
but only verbal explanation, literally: "because of the huge ab- 
sorption in Mg ii, the local continuum for the Mn i UV lines 
changes. There is less flux and thus fewer photons to be ab- 
sorbed and hence the ground level is consequently more popu- 
lated". What actually happens in such coupling is that the quasi- 
continuous Mgn wing opacity increases the height of photon 
escape and enforces LTE behavior to considerably larger height 
than the Mn i lines might maintain on their own, up to the height 
where Mg n k photon loss causes source function departure from 
the Planck function. Thus, thanks to the wavelength coincidence, 
departures from LTE set in only at exceptionally large height, 
and these anyhow represent a much larger fractional population 
change for the upper level than for the ground level, affecting the 
source function much more than the opacity. 



Optical pumping, also claimed in the paper's verbal ex- 
planation, is perhaps easier grasped. Super-Planckian radiation 
in a pumped transition may overpopulate its upper level and 
so induce super-Planckian source function excess and apparent 
brightening in subordinate lines from the same upper level, in 
this case Mn i multiplet 4, and perhaps also in other lines through 
upper-level coupling. Such pumping was earlier established for 
a var iety of emission lines in the extended w i ngs of Can H & K 
re.g.. lCanfieldlll971l: iRutten & Stengel 19801: ICram et al.lfT98Qh . 
but these arise through coupling to more deeply escaping super- 
Planckian radiation outside the H&K wings, a wholly differ- 
ent mechanism. No such pumping affects the Mm lines, no r 
did pumping operate in the computation of IDoyle et alJ (1200 ll) : 
the slight brightening of Mm 5394.7 A which they obtained for 
deep-lying chromospheres was simply contributed through the 
outer tail of the contribution function. 

Actually, the computation was intrinsically wrong because 
the Mn i multiplet UV1 lines do not coincide with the Mg n h & k 
cores but lie beyond AA = 0.7 A in their wings, of which the in- 
tensity is considerably overestimated when assuming complete 
redistribution instead of partially coherent scattering - one of 
the worst locations in the whole solar spectrum to make this 
mistake. In the computations deep onsets of the chromospheric 
temperature rise resulted in bright extended Mg n h & k wings, 
but in reality the independent radiation fields in the inner wings 
decouple from the Planc k function already in the photosphere 
dMilkev & Mihalas 1974). This error i s obvious when co mpar- 
ing the computed profiles in Fig. 3 of IDoyle et alJ (1200 ll) with 
observed Mgn h&k profiles such as the pioneering ones by 
iLemaire & Skumanichl(ll973l) in which even the strongest plage 
emission shows deep ki minima at AA = ±0.5 A from line cen- 
ter. In fact, the Mn 1 2794.82 A line is visible as an absorption dip 
within the deep Mgnki minimum not only in the reference spec- 
trum in Fig. 1 of IStaath & Lemaird (1 19951) but also in all Mgn 
h&k profiles displayed in their Fig. 2, and it remains located 
within the ki dip even in all limb spectra in their Fig. 11, both 
inside and outside the limb. The latter result from summation of 
chromospheric Mg ii k emission along the line of sight so that 
their peak widths represent a maximum. Therefore, everywhere 
across the solar surface Mm 2794.82 A lies in the deep ki dip 
which has a sub-Planckian source function due to coherent scat- 
tering and does not at all respond to chromospheric activity. 

In summary, although both the Mm 5394.7 A line and the 
peaks in Mgn h&k are observed to track solar activity, the 
blending of the Mn i UV1 lines into the opaque Mg n h & k wings 
do es not imply a viab le causal relationship, nor was one proven 
bv lDovle etafl (l200ll) . The Mm UV1 lines lie too far from the 
chromospheric h&k emission peaks for any pumping by these. 
In addition, the Mm multiplet 1 lines are formed much deeper 
(IVitas & Vincell2007l) . 

1.3. Photospheric interpretation 

Height of formation estimates for the optical Mn i lines based 
on standard modeling suggest s that they are purely photospheric 
([Gurtovenko & Kostvkl Il989l: IVitasI l2005h . Observational evi- 
dence tha t Mm 5394 . 7 A is p urely photospheric has been col- 
lected bv lVince et all (l2005bl) who showed that Mm 5394.7 A 
bisectors show charact eristic photospheric s hapes and center- 
limb behavior, and by iMalanushenko et all (120041) who com- 
pared spectroheliogram scans in Mm 5394.7 A and the nearby 
Mm 5420.4 A line with other lines and a magnetogram. The 
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Table 1. Line parameters. 



Fig. 1. Schematic G-band bright-point enhancement in a mag- 
netic concentration, portrayed here as a vertical cut through an 
idealized flaring fluxtube (thick curves) containing strong field, 
embedded in field-free quiet photosphere. Left: in radial view- 
ing the Wilson depression due to magnetic-pressure evacuation 
deepens the photon escape layer characterized as r = 1 surface 
to far below the outside surface. The magnetic concentrations 
are much cooler there than the subsurface surroundings due to 
suppressed convection, but hotter than the outside photon es- 
cape layer due to the large depth and hot-wall irradiation, and 
have flatter temperature gradients. The correspondingly larger 
degree of dissociation of the CH molecules that make up the G 
band cause a yet larger effective Wilson depression in this spec- 
tral feature, producing larger brightness enhancement than in the 
continuum. Right: In near-limb viewing the same lack of opac- 
ity along the slanted line of sight causes deeper, facular "bright- 
stalk" ^samrjlirmofhot granules behin d magnetic con centrations. 
From lRuttenl (119991) : see also Fig. 1 of ISpruitl (Il976b . 



Mn i lines show network bright, closely mimicking the unsigned 
magnetogram signal which is purely photospheric. 

The arguments above against a chromospheric interpretation 
and these observational indications of photospheric formation 
together suggest that the propensity of Mn i lines to track solar 
activity in their line-center brightness may be akin to the con- 
trast brightening that network and plage show in the G band. 
We therefore summarize G-band bright-point formation, where 
"bright points" stands for kilo-Gauss magnetic concentrations. 
Their enhanced photospheric brightness in continuum intensity 
and further contrast increase in G-band imaging has been stud- 
ied in extenso and is well understood, both for their on-disk ap- 
pearance as filigree and near-limb appearance as faculae. A brief 
review with the key references is given in the introduction of 
iDe Wijn et al.l(l2005l) : the carto on in Fig. [1] sch em atizes the mag - 
netostatic fluxtube paradigm o fHwaan (1967) and Spruitl d 1976b . 
The latest num erical verifications o f this concept a r e the MHD 
simulations of iKeller et all (120041) . ICarlsson et all (120041) . and 
IShelyag et al.1 (120041) . The upshot is that the G band shows mag- 
netic concentrations with enhanced contrast because its consid- 
erable line opacity lessens through molecular dissociation within 
the concentrations while its LTE formation implies good temper- 
ature mapping. Similarly, the extended blue wing of Ha bright- 
ens also in magnetic concentrations th rough lessening of coll i- 
sional broadening plus LTE formation dLeenaarts et al.ll2Q06bl) . 

Hence, for manganese we seek a property that enhances the 
reduction of line opacity in fluxtubes over that in comparable 
lines say from Fe i, enhancing the corresponding "line gap" phe- 
nomenon. A first consideration is that Mn 1 5394.7 A is relatively 
sensiti ve to temperature, as pointed out already by lElste & Te ske 
(119781) . because it originates from a neutral-metal ground state. 
Such lines do not suffer from the cancellation in LTE response to 
temperature increase that lines from excited levels have through 
masking of higher- temperature perturbations by radially moving 
the formation of the line outward where it samples lower tern- 



Line 


Mn i 


Fei 


Wavelength [A] 


5394.677 


5395.215 


Transition 


a 6 S 5/2 -z*P° 1/2 


z 5 G° 2 -g 5 F l 


Excitation energy [eV] 


0.0 


4.446 


Oscillator strength (log gf) 


-3.503 


-1.74 


Lande factor 


1.857 


0.500 


Slower 


-2.41 




^upper 


18.23 




Ionization energy [eV] 


7.44 


7.87 


Abundance (log Af H = 12) 


5.35 


7.51 



peratures (see Fig. 4 of iLeenaarts et al.ll2005|) . In addition, Mm 
5394.7 A is also a somewhat forbidden intersystem transition 
which makes its source function obey LTE more closely than for 
higher-probability lines. However, both these properties are un- 
likely to play an important role since Malanushenk o et a D (120041) 
found that Mn 1 5420.4 A, a member of multiplet 4 at 2. 14 eV ex- 
citation, brightens about as much in network. 

The obvious remaining property which makes Mn i lines dif- 
fer from others is their large hyperfine structure. How can this 
cause unusual brightness enhanceme nt in strong-field magnetic 
concentrations ? lElste & Teskel (Il978l) pointed out that it lessens 
sensitivity to the "turbulence" that was needed to explain other 
lines in classical one-dimensional modeling of the spatially- 
averaged solar spectrum. The ill-famous "microturbulence" and 
"macroturbulence" parameters were supposed to emulate the re- 
ality of convective and oscillatory inhomogeneities which make 
a solar-atlas line profile represent a spatio-temporal average over 
widely fluctuating and Dopplershifted instantaneous local pro- 
files. Lines that should be deep and narrow are so smeared into 
shallower average depressions. However, lines that are already 
wide intrinsically through hyperfine broadening suffer less shal- 
lowing by being less sensitive to Dopplershifts. The culprit may 
therefore not be the hyperfine structure of Mn i lines but rather 
the heavy hydrodynamic smearing of all the other, narrow lines 
that occurs in the granulation outside magnetic concentrations. 
We test this idea below and find it is correct, turning the analysis 
into one of general Fe i line formation rather than specific Mn i 
line formation. 

We first demonstrate this idea with classical one-dimensional 
modeling in Sect. 13.11 then verify it through three-dimensional 
MHD simulation in Sect. 13.21 and add explanation through dis- 
section of the simulation in Sect. 13.31 The next section presents 
our assumptions, input data, and methods. 



2. Assumptions and methods 

2.1. Line selection 

We started this project with a wider line selection but for the sake 
of clarity and conciseness we limit our analysis to Mn 1 5394.7 A 
and t he neighboring Fei 53 95.2 A line, following the exam- 
ple of iDanilovic et alJ (I2007I) who show and modeled the activ- 
ity modulation of both lines in parallel. The Fei 5395.2 A line 
serves here as prototype for all comparable narrow lines. The 
line parameters are given in Table [T] The Mn i hyperfine struc - 
ture constants J{ come from Blackwell- Whitehead et al. d2005h . 
All other values were taken from the NIST database at URL 
http : //physics . nist . gov/ asd3, 
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Fig. 2. One-dimensional standard models plotted as tempera- 
ture against column mass. The lower soli d curve is the quiet-sun 
MACKKL model of iMaltbv etaD (Il986l) : the upper solid curve 
is the PL A model for fluxtubes in plage of ISolanki & Brigljevid 
(I1992I) . The tick marks on these two curves specify = 1 loca- 
tions for the line-center wavelength of Mm 5394.7 A (upward 
arrows), the line-center wavelength of Fe i 5395.2 A (downward 
arrows), and the continuum in between (arrowless). Extra Tx - 1 
ticks are added for line synthesis without turbulent smearing, but 
these differ significantly only for the F ei line and MACKK L. 
The two dotted models were defined by lUnruh et all (1 19991) as 
basis for solar irradiance modeling and are discussed in Sect. |U 
The lower one, close to MACKKL, is for quiet sun, the upper 
one for plage. 



2.2. Line synthesis 

We perform line synthesis in the presen c e of m agnetic fields 
with the code of Sanchez Almeida et al. I <l2008h which solves 
the radiative transfer equation for polarized light in a given 
one-dimensional atmosphere via a predictor-corrector method. 
It yields the full Stokes vector, but we only use the inten- 
sity here. The code includes evaluation of the Zeeman pat- 
tern for lines with hyperf ine structure using the routine of 
lLandi degrinnocentil (Il978l) . This pattern depends on the mag- 
netic field and on the hyperfine structure constants, the quantum 
numbers of the upper and lower level, the relative isotopic abun- 
dance, and the isotope shifts. The splitting depends on the hy- 
perfine structure constants {A and £>, which account for the two 
first terms of the Hamiltonian describing the interaction between 
the electrons in an atomic level and the nuclear magnetic mo- 
ment. {A quantifies the magnetic-dipole coupling, & the electric- 
quadru pole coupling. We consider B n egligible here, again fol- 
lowing ls"anchez Almeida et aD d2008) who successfully repro- 
duced multiple Mn i line profiles with different HFS patterns. 

2.3. Assumption of LTE 

The source function of Fe 1 5395.2 A is likely to share the charac- 
teristic properties of weak subordinate Fe i lines to possess LTE 
source functions (with the equilibrium maintained by the much 
stronger Fe i resonance lines) while having less-than-LTE opac- 
ity in the upper photosphere in locations with steep temperature 




0.1 r ~ 

0.0 Ej , , , i , , , i , , , i , , , i , , , i , c 

5394.4 5394.6 5394.8 5395.0 5395.2 5395.4 
A [A] 

Fig. 3. Upper panel: line extinction coefficient a\ for Mm 
5394.7 A (left) and Fei 5395.2 A (right), for temperature T = 
5000 K and longitudinal magnetic field strength B = G (top 
solid curve), and 2000 G (lowest solid curve). For the Mm line 
there is also an intermediate profile for B = 1000 G. The dotted 
curves result when temperature T = 7000 K is inserted into the 
Dopplerwidth. The insets specify the splitting pattern of Mn i 
5394.7 A for B = G (upper) and B = 2000 G (lower). Lower 
panel: line profiles resulting from classical one-dimensional 
spectral synthesis. The lower solid curve in th e lower trio is 
the quiet- sun disk-center spectrum in the atlas of lNeckell (1 19991) 
which is an absolute -intensity version of the NSO/FTS atlas of 
I Wallace et al.l (Il998l) . The other solid curve in the lower trio 
is the spectrum computed from the MACKKL model includ- 
ing micro- and macroturbulence. The dotted curve results when 
these fudge parameters are omitted. The two upper spectra result 
from the PL A model with B = 1000 G, also with and without 
turbulent smearing. 



gradients (as controlled by rad iative oyerion ization of iron in the 
near ultraviolet: see review byEutten 198j). 

The source function of Mm 5394.7 A should also remain 
rather close to LTE because it is an intersystem transition 
with rather small oscillator strength, although not as forbid- 
den as the well-known LTE Mgi 4571.1 A line. Indeed, Mm 
53 94.7 A closely obeys LTE in t he detailed NLTE computations 
by lBergemann & Gehrenl d2007) for a standard solar atmosphere 
model assuming radiative equilibrium. Its opacity is set by the 
ground-state population which is close to LTE everywhere in 
these computations, even when photoionization is increased by a 
factor 5000. The degree of manganese ionization is likely to ex- 
ceed LTE in locations with steeper radial temperature gradients, 
but so will the degree of of iron ionization; there is no reason 
to suspect large difference between Mn i and Fe i line formation 
with respect to NLTE effects. Thus, apart from their large hyper- 



N. Vitas et aL: Explanation of the activity sensitivity of Mm 5394.7 A 



5 



fine structure, Mn i lines should not behave very differently from 
Fe i lines. 

The radial temperature gradients within magnetic concentra- 
tions are probably c lose to local radiative eq uilibrium throughout 
their photospheres ( Shemmovaetal 1120051) . so that NLTE overi- 
onization is likely to affect both lines similarly also in these. 

2.4. One-dimensional modeling 

We use the standard model PLA for fluxtubes making up 
plage that was derived in the 1980s by Solanki and cowork- 
ers f r om spectropolarimetrv of photospheric lines dSolankil 
11986b ISolanki & Steenbockl Il988l ISolanki & Brigljevicl Il992h 



and shown in Fig. 1 of lBruls & Solankil(ll993l) . Since we only in- 
tend demonstration here, we do not apply spatial averaging over 
upw ard-expanding and c anopy-merging magnetostatic fluxtubes 
as in lBunte et al.1 (Il993l) . but simply use it as on-axis represen- 
tation of a fully-resolved fluxtube as in the cartoon in Fig.[T] 
with consta nt field strength along th e axis. The same was done 
in Fig. 4 of Sanchez Almei da et al. PLA is shown in 

Fig. [2] togethe r with the standard MACKKL quiet- sun model of 
Malt bv etal.1 (119861) which we use to represent the non-magnetic 
atmosphere outside fluxtubes. In a plot like this PLA appears 
to be much hotter than MACKKL, but when PLA is shifted 
over the Wilson depression in fluxtube modeling it is actually 
much cooler at equal geometrical height. Figure |2] shows two 
additional models from lUnruh et al.l ([1999) that are discussed in 
Sect. [4] 

2.5. Three-dimensional simulation 

We use a single snapshot from a time-dependent simulation with 
the MURaM (MPS/Univers i ty of Chicago Radiative MHD) cod e 
(IVogler & Schusslerl 120031: IVoglerl 120041: IVogler et al.l 12005b . 
MURaM solves the three-dimensional time-dependent MHD 
equations including non-local and non-grey radiative transfer 
and accounting for partial ionization. 

The particular snap shot used here is tak en from a simulation 
as the one described by|V6 gler et al. I d200l) . Its horizontal extent 
is 6 Mm sampled in 288 grid points per axis, its vertical extent 
1.4 Mm with 100 grid points. It was started with a homogeneous 
vertical seed field of 200 G. Magnetoconvection gave it an ap- 
pearance similar to active network with strong-field magnetic 
concentrations in intergranular lanes. A relatively quiet subcube 
was selected for our line synthesis. It has 2.5x2.5 Mm horizontal 
extent and contains a large granule and a field-free intergranular 
lane, in addition to lanes containing magnetic concentrations of 
varying field strength (Fig. [4]). 

The line synthesis performed for this paper used the code 
described above, treating the vertical columns in the subcube as 
independent lines of sight. The displays below are restricted to 
the nominal NIST wavelengths of the two lines in Table [T] cor- 
responding to the line centers in a spatially-averaged disk-center 
intensity atlas as in Fig. [3] 

3. Results 

3.1. One-dimensional demonstration 

Figure [3] presents results of our one-dimensional modeling. The 
solid curves in the upper panel show the spectral variation of the 
extinction coefficient for both lines for temperature T = 5000 K 
and field strengths B = G and B = 2000 G, plus the inter- 
mediate profile for B = 1000 for the Mn i line. For the Fe i line 



the Zeeman effect produces simple broadening but for the Mn i 
line the many hyperfine components, each with its own magnetic 
splitting, cause an intricate pattern shown in the lower part of the 
inset. The resulting profile widens in the wings but the core first 
becomes peaked at B = 1000 G and then splits into three collec- 
tive peaks at B = 2000 G. 

The dotted curves result from inserting temperature T = 
7000 K into the Dopplerwidth but not into other variables, in or- 
der to show the effect of larger thermal broadening while keep- 
ing all other things equal. The Fe i profile for B = G looses 
appreciable amplitude while the Mn i profile does not. 

The lower panel of Fig. [3] shows emergent intensity pro- 
files for the two lines. The lowest solid curve is the observed 
spatially-averaged disk-center spectrum. It is closely matched 
by the MACKKL modeling when applying standard microtur- 
bulence (1 km s" 1 ) and best-fit macroturbulence (1.28 km s" 1 
for Mm 5394.7 A, 1.55 km s" 1 for Fei 5395.2 A). When this ar- 
tificial broadening is not applied the computed Fe i line becomes 
too deep, but the depth of the Mn i line does not change thanks 
to its flat-bottomed core. The upper curves result from the PLA 
model with B = 1000 G, again with (solid) and without (dotted) 
turbulent smearing. The smearing again affects only the Fe i line. 
It causes a corresponding shift of the Fe i = 1 location along 
MACKKL in FigS 

Comparison of these MACKKL and PLA results shows that 
the spectrum brightens at all wavelengths but most in the Mn i 
line, by a factor 2. Turbulent smearing does not affect this line 
but it produces large difference for the Fei line. In particular, 
if it is applied to the MACKKL quiet- sun prediction but not to 
the PLA profile, the Fe i line-center brightness increase is only a 
factor 1.4. 

This difference in line-center brightening suggests that the 
Fei line suffers more from thermal broadening and the ther- 
modynamic fine structuring that was traditionally modeled with 
micro- and macroturbulence. The apparent sensitivity of the Mn i 
line to magnetic activity may therefore indeed result from non- 
magnetic quiet-sun Doppler smearing of the Fe i line, lessening 
the latter' s sensitivity. 



3.2. Three-dimensional verification 

The results from the MURaM simulation are shown in Fig.SJ 
The three panels in the upper row show basic state parameters 
across the T5 = 1 surface where T5 is the continuum optical depth 
at A = 5000 A determined separately for each simulation column. 
The middle row displays synthetic intensity images for our three 
diagnostics: the continuum between the two lines and the nomi- 
nal line-center wavelengths of Fe 1 5395.2 A and Mn 1 5394.7 A. 
The bottom row shows these intensities in the form of pixel-by- 
pixel scatter plots against the magnetic field strength at T5 = 1 . 

The three images demonstrate directly why Mm 5394.7 A 
shows larger brightness contrast between non-magnetic and 
magnetic areas: the magnetic bright points reach similar bright- 
ness in all three panels but the granulation is markedly darker in 
this line. Darker granulation implies larger sensitivity to activity, 
i.e., addition of more magnetic bright points. 

The three scatter plots in the bottom row of Fig. 01 quantify 
this behavior. In the continuum plot at left the darkest pixels lie 
in field-free or weak-field intergranular lanes. Addition of mag- 
netic field within the lanes brightens them, more for stronger 
fields, and for the strongest fields almost up to the maximum 
brightness of field-free granular centers. 
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Note that in slanted near-limb viewing the magnetic con- 
centrations do not add brightness to dark intergranular lanes but 
permit deeper viewing into bright granules behind them, adding 
brightness to the already brightest features and so making facu- 
lae brighter than the granular background (see Fig. [I}. 

The scatter diagram for Fe i 5395.2 A (bottom-center panel) 
shows a similar hook shape. The cloud of granulation pixels at 
left lies lower since the line is an absorption line. However, the 
pixels with the strongest field still brighten to the same values 
as in the continuum panel, implying that the line vanishes com- 
pletely at its nominal wavelength. 

The scatter diagram for Mn i 5394.7 A shows similar behav- 
ior, but it loses the hook shape. The line is yet darker in field-free 
granulation. In this case the corresponding dark cloud of low- 
field pixels at left does not have an upward tail. However, the 
upward tail of pixels with increasing field still stretches all the 
way from the dark lanes to the continuum values, which again 
implies line vanishing. 

If magnetic field is added to field-free granulation, this ad- 
dition takes away points from the bottom left of the distribution 
(corresponding to dark intergranular lanes) and adds "magnetic 
bright points" at the upper right. Its effect is more dramatic in 
the ensemble average of the Mn i line than in that of the Fe i line 
because the latter already has bright contributions from field- 



Fig. 4. Results from the MURaM snapshot. 
First row: magnetic field strength, temperature, 
and gas density across the r 5 = 1 continuum- 
escape surface. Greyscale calibration bars are 
shown on top. The field is monopolar. The 
white line specifies the cut used in Figs. \5\- 
[8] to display formation parameters in a ver- 
tical plane through the simulation volume. It 
represents a spectrograph slit in Fig.|7] The 
three superimposed ticks specify locations an- 
alyzed in detail in Fig. 13 a magnetic concen- 
tration, an intergranular lane, and the edge of 
the large granule in the center. Second row: 
synthetic intensity images for the continuum 
between the two lines and the nominal wave- 
lengths of Fei 5395.2 A and Mm 5394.7 A. 
They share the same intensity greyscale. These 
images represent predicted observations with a 
telescope matching the MURaM resolution by 
being diffraction-limited at 2.4-m aperture. The 
magnetic concentrations (point-like locations 
with the largest field strength) are about equally 
bright at all three wavelengths, but the non- 
magnetic granules (in the center and at the top 
of the field) are much darker in Mm 5394.7 A 
than in the other two panels. Third row: the 
same intensities plotted as scatter diagrams per 
pixel against the magnetic field strength at the 
r 5 = 1 level for that pixel. The upward tails of 
bright values at large field strength reach the 
same maxima in all three plots, but the gran- 
ulation (B < 200 G) is much darker in Mm 
5394.7 A than in Fei 5395.2 A and the contin- 
uum. The first two panels display a pronounced 
hook shape, but this hook is not present in the 
Mn i panel. Adding more field therefore gives 
more spatially- averaged brightening in the Mn i 
line than in the Fe i line and the continuum. The 
Mm line thanks its "unhooking" to its hyper- 
fine structure through which it is less sensitive 
to the Doppler brightening that affects the Fe i 
line (and all similar lines) in the granulation. 

free granules; in the Fe i line the field-free granulation covers the 
same brightness range as the field-filled lanes. Thus, field addi- 
tion means larger spatially-averaged brightness increase in the 
Mn i line than in the Fe i line. 

We conclude that the MURaM synthesis duplicates the sun 
in having larger brightness response to more activity in the Mn i 
line than in the Fe i line - not via chromospheric emission but 
through showing granulation darker. 

3.3. Explanation 

Unlike the sun, the MURaM simulation offers the opportunity 
to not only inspect the emergent spectrum but also to dissect the 
behavior of pertinent physical parameters throughout the simula- 
tion volume. So instead of ending this paper here with the above 
conclusion that neither a chromosphere nor NLTE coupling to 
Mg ii h & k are needed to reproduce the activity sensitivity of 
Mm 5394.7 A, we add analysis of the MURaM results to diag- 
nose why the granulation appears darker in the Mm line - or 
rather, why it appears brighter in the Fe i line - and so earn our 
use of the term "explanation" in the title of this paper. 

The white line with three ticks in the grey- scale panels of 
Fig. H] specifies the spatial samples that are selected in Figs.[5]-[8j 
This particular cut and the tick locations were chosen to sam- 
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Fig. 5. Various quantities in the MURaM snapshot across the 
vertical plane marked by the white horizontal line in Fig.|U The 
upper panels show the magnetic field strength and radial veloc- 
ity as a function of geometrical height versus horizontal location 
along the cut. The height scale has h = km at the mean T5 = 1 
level, averaged over all columns. It extends 882 km below and 
518 km above that level in the simulation but only the perti- 
nent height range is shown here. The velocity coloring is red for 
downdraft (redshift) and blue for updraft (blueshift). The lower 
panels show the temperature and the gas density. The overlaid 
curves specify T5 = 1 continuum formation depths. The three 
ticks mark the locations analyzed in Fig. [7] 



pie a strong-field magnetic concentration that appears as a bright 
point in the intensity images (lefthand tick), a non-magnetic dark 
lane (middle tick), and a granule (righthand tick). Unfortunately, 
this cut samples only the edge rather than the center of the large 
granule covering the center part of the field but otherwise we 
wouldn't have sampled both a bright point and a non-magnetic 
lane. A magnetic lane less extreme than the bright point is sam- 
pled at the far right. 

Figure [5]diagnoses MURaM physics in the vertical plane de- 
fined by this cut. It shows behavior that is characteristic of solar 
magnetoconvection near the surface. The overlaid curves spec- 
ify the T5 = 1 continuum surface. The very low gas density (4th 
panel) in the magnetic concentration at x « 400 km produces a 
large Wilson depression of about 200 km. The deep dip in the 
T5 = 1 curve so samples relatively high temperature, as well as 
a relatively flat vertical temperature gradient. The intergranular 
lane at x « 1300 km combines low temperature with high den- 
sity and strong downdraft; the granule edge at x « 1650 km 
combines higher temperature with gentler updraft and low sub- 
surface density. 

Figure [6] repeats this vertical-plane display of temperature 
and density but plotted per column on the radial optical depth 
scale that belongs to each diagnostic. The panels in the first and 
second columns show the atmosphere "as seen" by each spectral 
feature at its nominal wavelength. The dotted horizontal lines at 
logr^ = indicate their formation heights. The solid curves are 
the T5 = 1 locations. The third column shows relative behavior of 



the corresponding opacities in the form of fractional lower-level 
population variations. 

The top panels of Fig.[6]show only slight differences between 
the r = 1 locations along the cut. The magnetic concentration has 
an appreciable hump in T(tx) and dip in p(jx) around T5 = 1. 
The relatively high temperature and low density there combine 
into increased electron-donor ionization. This is illustrated by 
the sixth panel showing the fractional population variation of 
the lower level of Mn 1 5394.7 A, which is the Mn 1 ground state. 
Its behavior equals the depletion by manganese ionization apart 
from a minor correction for the temperature sensitivity of the 
Mn 1 partition function. Manganese has too small abundance to 
be an important electron donor, but it ionizes similarly to iron 
(Table [B so that this panel illustrates characteristic electron- 
donor ionization, with neutral-stage depletion occurring within 
magnetic concentrations and at large depth. The corresponding 
increase of the free-electron density produces larger H" opac- 
ity, evident as overall color gradient reversal between the top 
and center panels in the third column. It results in an upward 
enhancement peak at the magnetic concentration in the FT pop- 
ulation panel. This relative increase of the continuum opacity 
explains that Tx = 1 is reached at lower density in the mag- 
netic concentration than in the adjacent intergranular lane (sec- 
ond panel); the Wilson depression is smaller than one would es- 
timate from pressure balancing alone. The flat temperature gra- 
dient in the magnetic concentration produces a marked upward 
extension in T(tx) in the first panel. It contributes brightness en- 
hancement along much of the intensity contribution function and 
so makes magnetic concentrations appear bright with respect to 
non-magnetic lanes. However, the hook pattern in the continuum 
scatter plot in Fig. [4] shows that this lane brightening does not ex- 
ceed granular brightness. 

The center row of Fig. [6] shows that Mm 5394.7 A is gen- 
erally formed higher than the continuum, but not in the mag- 
netic concentration where the upward hump in the T5 = 1 curve 
nearly reaches the 7^ = 1 level. This is because the relatively 
high temperature and low density there increase the degree of 
manganese ionization, as evident in the third panel which shows 
a marked upward blue extension. The line weakens so much that 
the brightest Mm pixels in Fig. |4] reach the same intensity as in 
the continuum, also sampling the upward high-temperature ex- 
tension (first column). Conversely, the largest line-to-continuous 
opacity ratio (i.e., the largest separation between the = 1 and 
T5 = 1 locations) is reached along the intergranular lane where 
the neutral- stage population is larger due to relatively large den- 
sity and low temperature. Thus, the curve separation maps the 
variation of the fractional ionization along the dotted line. Since 
the temperature increases inward anywhere, the deeper sampling 
within the magnetic concentration and the higher sampling in the 
granulation, especially in the lanes, together increase the bright- 
ness contrast between bright point and granulation compared to 
that in the continuum. This effect of increased ionization is sim- 
ilar to the effect of CH dissociation in the G band and reduced 
damping in the Ha wings within magnetic concentrations. It en- 
hances the magnetic lane brightening so that that exceeds granu- 
lar brightnesses, undoing the hook shape of the continuum scat- 
ter plot. 

The bottom row of Fig. [6] shows the corresponding plots for 
Fei 5395.2 A. Iron and manganese ionize similarly so that low- 
excitation Fe 1 lines suffer the same depletion in magnetic con- 
centrations. However, Fei 5395.2 A is a high-excitation line; 
its Boltzmann sensitivity to higher temperature largely compen- 
sates for the enhanced ionization so that the panel in the third 
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Fig. 6. Temperature, gas density, and popula- 
tion variations across the vertical cut through 
the MURaM simulation plotted per column on 
radial optical depth scales for the nominal con- 
tinuum (top row), Mn 1 5394.7 A (middle row), 
and Fei 5395.2 A (bottom row) wavelengths. 
The solid curves specify r 5 = 1 depths, the 
horizontal dotted lines feature- specific r A = 1 
depths. The third column displays fractional 
population offsets. These are the lower-level 
populations for FT bound-free transitions, the 
Mm line, and the Fei line, each normalized 
by the total element density (hydrogen, man- 
ganese, iron) per location and shown in loga- 
rithmic units scaled to the mean value at t a = 1 
along the cut. Positive offsets are colored am- 
ber, negative bluish. The scale runs from -1.8 
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column looks much more like the top one for H" than like the 
second one for Mm 5394.7 A: the overall color gradient flips 
back again. The line again vanishes in the magnetic concentra- 
tion, so that the strongest- field pixels in Fig.|4] again reach the 
continuum intensities. One might expect that the curve sepa- 
ration in the bottom panels would follow the temperature pat- 
tern under the dotted line in the first column, but this is not 
the case; for example, the line also nearly vanishes (at its nom- 
inal wavelength) in the intergranular lane. This discordant vari- 
ation is caused by the Dopplershifts imposed on the line extinc- 
tion by the flows displayed in the upper-right panel of Fig.[5j 
Comparison shows that the unsigned amplitude of the flow vari- 
ation along the cut at 100 - 200 km above the T5 = 1 curve 
there is mapped precisely into reversed modulation of the T5 = 1 
curve for Fei 5395.2 A here. Thus, the line vanishes because it 
is shifted aside from its nominal wavelength. 

We demonstrate the Doppler-related formation differences 
between the two lines further in Fig. [TJ The upper two panels 
display spectral representations along the cut defined in Fig.Hl 
in the form of a spectral image and of profile samples. The top 
panel shows what would have appeared in an observational spec- 
trogram from a telescope with MURaM resolution. Both lines 
show a bright "line gap" in the magnetic concentration near the 
bottom. At x « 2500 km the intergranular lane with less strong 
field still causes a noticeable gap. Large Dopplershifts occur in 
the field-free lane and the granular edge. 

The second panel shows the product of data reduction of 
this spectrogram. Both panels illustrate that the line-center in- 
tensity of the Mm line is not very sensitive to Dopplershifts, 



whereas the Fei line shifts well away from its nominal wave- 
length nearly everywhere. It also weakens more in the granule 
edge from larger thermal broadening (Fig. [3]). Taking the spatial 
mean at each nominal wavelength does a fair job of intensity av- 
eraging for the Mn 1 line but misses nearly all dark cores in the 
Fei line, especially in the lanes but also in granules. Note that 
this particular cut does not represent hot granules well; more 
samples of these would add many profiles with weakened cores 
blueward of the nominal Fe 1 wavelength. 

The third panel shows the spectral profiles for the three sam- 
ple locations along the cut. The nominal Fe 1 line-center wave- 
length misses all three cores! Thus, the brightness average at this 
wavelength is much higher than it would be for an undisturbed 
line of this opacity. The Mn 1 line-center wavelength, however, 
only misses the deepest part of the magnetic-concentration pro- 
file which is weak anyhow. This disparity in Doppler sensitivity 
explains why the granulation in Fig.|4]is much darker in Mm 
5394.7 A than in Fei 5395.2 A. 

Also note the sharpening of the Mn 1 line profile from boxy 
to more pointed in some of the intermediate profiles in the sec- 
ond panel, which follows the extinction coefficient behavior in 
Fig.[3l It contributes to the decrease in equivalent width of Mm 
5394.7 A at larger activity, which is not analyzed here. 

The bottom panels of Fig. [7] show the vertical temperature 
stratifications at the three sample locations, at left against geo- 
metrical height with h = at T5 = 1 for the simulation mean, 
at right against column mass per feature. These graphs link the 
simulation results back to the classical fluxtube modeling in 
Sect. 13.11 They display familiar properties of granulation and 
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Fig. 7. Further analysis of the formation of 
Mm 5394.7 A, Fei 5395.2 A, and the inter- 
mediate continuum in the MURaM simulation. 
Top panel: spectrum synthesized from the com- 
puted emergent intensities. The vertical coordi- 
nate corresponds to the x sampling along the 
cut specified in Fig.|4] The three black hori- 
zontal markers correspond to the ticks select- 
ing the bright point (lowest), field-free inter- 
granular lane (middle), and granule edge (high- 
est). The three vertical lines specify the nom- 
inal wavelengths for Mn i 5394.7 A (left), Fe i 
5395.2 A (right), and the intermediate contin- 
uum (center). Second panel: spectral profiles 
for every second pixel along the cut. Third 
panel: spectral profiles for the three locations 
specified by ticks in Fig.[4] and black markers 
in the top panel. Dotted: magnetic concentra- 
tion. Dot-dashed: intergranular lane. Dashed: 
granule edge. Bottom panels: temperature strat- 
ifications for the three selected locations, at left 
against geometrical height, at right against col- 
umn mass as in Fig. [2 The solid curve is the 
PLA model, as in Fig. [2 The tick marks spec- 
ify the Tx — 1 locations for the nominal wave- 
lengths of Mm 5394.7 A (upward arrows), of 
Fei 5395.2 A (downward arrows), and of the 
continuum in between (arrowless). In the mag- 
netic concentration the Fei 5395.2 A tick co- 
incides nearly with the continuum one. The 
Wilson depression is much larger for the Mn i 
line than for the continuum and the Fe i line. 



magnetic concentrations: the temperature gradient is steepest in 
granules, flatter in intergranular lanes, with a mid-photosphere 
cross-overproducing "reversed granulation", and flattest in flux- 
tubes. The heights of formation, marked by ticks specifying the 
r A = 1 locations, are similar per line in the granular edge and 
the lane but much deeper in the magnetic concentration which 
is the coolest feature at equal geometrical height but the hottest 
at equal column mass and at the T con t = 1 location per feature. 
The Wilson depression is about 200 km in the continuum and 
twice as large for the Mm line through enhanced ionization 
within the magnetic concentration, but not for the Fe i line due to 
its Doppler-deepened granulation sampling. Both lines weaken 
severely in the magnetic concentration through this deep sam- 
pling, and both are redshifted away from their nominal wave- 
length, Fe i 5395.2 A all the way so that its Tx = 1 tick coincides 
with the continuum one. 



4. Discussion 

4.1. Comparison of 1D modeling and 3D simulation 

We have used both classical empirical fluxtube modeling 
(Sect. l3.lt and modern MHD simulation (Sect. l3.2l) . The bottom 
panels of Fig. [7] connect the latter with the former. The right- 
hand panel shows remarkably close agreement between the PLA 
model and the MURaM magnetic-concentration stratification. 
This may be taken as mutual vindication of these very different 
techniques, but not as proof of correctness since both assumed 
LTE ionization (discussed below) while the sun does not. 

The simulation resolves the granular Dopplershifts that 
were emulated by turbulence in the classical modeling. The 



simulation analysis in Sect. 1 3. 3 1 confirms that Mm 5394.7 A 
activity-brightens more than Fei 5395.2 A because its intrinsic 
hyperfine- structure broadening suppresses thermal and Doppler 
brightening in non-magnetic granulation. This was already indi- 
cated by the ID demonstration in Fig. [3] 



4.2. Validity of LTE 



We have assumed LTE throughout this paper. In Sect. I2.3I we 
noted that the most likely departure from LTE which may affect 
our two lines is loss of opacity through radiative overionization, 
in locations with steep temperature gradients where the angle- 
averaged intensity exceeds the Planck functi on in ultraviolet 
photo-ionization edges (see e.g.. iRuttenl I2003I for more expla- 
nation). The steepest temperature gradients in Fig.[5] occur not 
along radial columns but transverse through the magnetic con- 
centration. Therefore, the ionization which contributes to the 
vanishing of the Mn i line from the corresponding bright point 
may be underestimated. If so, the H" opacity is equally underes- 
timated. However, the MURaM simulation itself assumes LTE 
with coarse spectral sampling and is likely to differ intrinsically 
if NLTE overionization affecting the FT opacity is taken into ac- 
count. Such overionization was similarly neglected in the em- 
pirical construction of the PLA model based on Fe i-line spec- 
tropolarimetry assuming LTE to define the optical depth scales. 
Improvement requires non-trivial evaluation of the ultraviolet 
line haze. 
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Fig. 8. The effect of excitation energy and line 
strength on Fei line formation. Each panel is 
the fractional population of the lower level of 
an artificial line comparable to Fei 5395.2 A 
along the cut through the MURaM simula- 
tion defined in Fig.Hl again shown on nomi- 
nal line-center optical depth scales, normalized 
to the mean value along the dotted line, and 
with the same logarithmic color coding as in 
the third column of Fig.[6j clipped at offsets 
-1.8 and +1.8. The solid curves are again the 
r 5 = 1 locations for the continuum at 5000 A, 
the dashed lines the line-center r A = 1 for- 
mation indication. Top row: the Fei 5395.2 A 
line at excitation energy 2, 3, and 4 eV, re- 
spectively. The corresponding mean Boltzmann 
factors at r A = 1 are 1.2 x 10" 3 , 1.7 x 10~ 4 , 
and 5.2 X 10" 5 . Second row: the Fei 5395.2 A 
at excitation energy 0, 3, and 6 eV but with 
the oscillator strength scaled to maintain the 
same line strength in the emergent spectrum 
that Fei 5395.2 A shows (by 25 x 10" 3 , 25 and 
25 x 10 3 , respectively). Bottom row: idem, but 
the color coding now measures deviations from 
the Milne-Eddington approximation by show- 
ing the lower-level population divided by the 
FT population. The leftmost panel is severely 
clipped. The Milne-Eddington approximation 
improves dramatically with higher excitaton. 



4.3. Extension to other lines 

We have demonstrated and analyzed Doppler brightening in de- 
tail only for Mm 5394.7 A and Fei 5395.2 A. The first may be 
taken to exemplify all lines with wide boxy extinction profiles, 
the second all lines with narrow peaked profiles, i.e., almost all 
other photospheric lines. Figure [8] serves to demonstrate the lat- 
ter generalization by extending the MURaM synthesis to a range 
of artificial Fe i lines by simply placing Fe 1 5395.2 A at other ex- 
citation energies. The first two rows again display relative pop- 
ulation variations across the simulation cut using the same color 
coding as in the third column of Fig. [6] The top row of Fig. [8] 
combines the effects of the Boltzmann factor on the amplitude 
of the line extinction and on its temperature sensitivity. The sec- 
ond row isolates the latter by making each line just as weak as 
Fei 5395.2 A. 

The first panel of Fig. [8] is included only for illustration be- 
cause this artificial line is already too strong for reliable synthe- 
sis from our simulation and would also suffe r appreciable de- 
partur es from LTE if it existed in reality (see iRutten & Kostikl 
Il982h . It is formed near the top of the simulation volume which 
sets the columnar structure in the upper part; the sizable vertical 
flows there impose the separation between the two r = 1 curves. 
In reality, very strong Fe i lines suffer less convective and thermal 
D oppler weakening in their cor es which explains the observation 
by lMalanushenko et al.l (12004b that the strongest Fe i lines show 
some (but only slight) global activity modulation. Even that 
modulation has nothing to do with the chromosphere! The sec- 
ond and third panel describe deeper line formation with steeper 



gradients. Their r = 1 separations also demonstrate Doppler 
modulation. Thus, Doppler weakening affects these three lines 
just as it affects Fe i 5395.2 A. 

The second row of Fig. [8] illustrates the effect of Boltzmann 
temperature sensitivity at given emergent line strength. The first 
panel describes the Fe i ground state and is nearly the same as 
the sixth panel of Fig. [6] Towards higher excitation energy the 
Boltzmann increase compensates more of the depletion by ion- 
ization. In the center panel the population offset gradients are 
very flat; at right, they are reversed so that the population varia- 
tion mimics the behavior of FT in the third panel of Fig.[6]closely, 
although at slightly smaller amplitude. The Doppler modulation 
of the r = 1 separation remains similar, irrespective of this gradi- 
ent reversal in fractional population distribution. Thus, Doppler 
brightening affects all photospheric Fe i lines similarly. 

The third row of Fig.[8]repeats this test in terms of the Milne- 
Eddington approximation (constant line-to-continuum opacity 
ratio with height ) which is often assumed in photospheric 
polarimetrv (e.p.. |Auer et all 119771: Iskumanich & Litesl 119871: 



IWestendorp Plaza et al.l 119981: lOrozco Suarez et al.l l2007b ~Al 
left the combination of Fe i depletion by ionization and the cor- 
responding increase in FT opacity implies squaring of the popu- 
lation offsets or twice steeper gradients in this logarithmic plot 
than in the second row. The Milne-Eddington approximation 
fails badly but least within the magnetic concentration. It im- 
proves with excitation energy through the Boltzmann compen- 
sation. The panel at right shows the residue from subtracting the 
sixth panel in this figure from the third panel in Fig.[6j resulting 
in rather small deviations. This behavior was shown earlier in 
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Fig. 3 of Rut ten & van der Zalm (1984; reprinted in Fig. 9.2 of 
lRuttenll2003h . 

4.4. Irradiance modeling 

In this paper we have not made the step from radial view- 
ing to full-disk averaging. The demonstration in Sect. 13. II used 
Solanki's one-dimensional empirical model for a fluxtube in 
plage but without adding flaring fluxtube geometry, granular 
presence next to a fluxtube and behind it in slanted facular view- 
ing, multiple-tube interface geometry, spatial averaging over 
these geometries, the multi-angle evaluation needed to emu- 
late center-to-limb viewing, and without full-disk averaging and 
consideration of spatial distributions over the solar surface and 
their variations with the solar cycle. All these non-trivial aspects 
would need careful quantification to expand the demonstration 
of enhanced sensitivity in Sect. 13.11 into a quantitative estimate 
for comparison with Livingston's irradiance data. The simula- 
tions in Sect. 13.21 yielded profile synthesis for a small area of 
solar surface containing strong-field concentrations that may be 
considered more realistic than idealized magnetostatic fluxtubes 
but, nevertheless, generation of full-disk signals comparable to 
Livingston's data still requires all the above quantifications. This 
effort is not done here. 

In contrast, the step from one-dimensional line synthesis to 
emulation of the full- disk and cycle-depende nt integrated signal 
was recently made by IDanilovic et al.1 (120071) using the SATIRE 
(Spectral A n d Tot al Irradiance R e constr uction) approach of 
Flig ge et alJ (l200Qh . iKrivova et all (120031) and Wenzler et al. 
(2005, 2006). In this technique the spatial distributions of spots 
and plage are extracted from full-disk magnetograms and contin- 
uum images to derive disk-coverage distributions throughout the 
more recent activity cycles. Each component (quiet sun, plage, 
spots) is represented by a standard one-dimensional model at- 
mosphere. The first two are shown as dotted curves in Fig. [2] 
The lower one for quiet sun is the radiative equilibrium model 
of Kurucz (1979, 1992a, 1992b) which is nearly identical to 
MACK KL. The upper one for plag e was made by lUnruh et al.l 
(119991) by smoothing model P of iFontenla et al.l (1 1993b and 
deleting its chromospheric temperature rise. For spots a Kurucz 
radiative-equilibrium model with low effecti ve temperature is 
used, not shown here. IDanilovic et al.l (120071) found that using 
these models with the empirically established SATIRE coverage 
fractions gives a good reproduction of Livingston's data both for 
Mm 5394.7 A and Fei 5395.2 A. 

Such use of the dotted models in Fig. [2 is an ad-hoc trick 
to reproduce the larger brightness of plage. The Mm 5394.7 A 
line activity-brightens more than Fei 5395.2 A in this simplis- 
tic modeling because it is a stronger line, hence formed higher, 
hence getting more out of the divergence between the two mod- 
els with height. Any line as strong would show the same bright- 
ening. Stronger lines would brighten more; in particular, Fei 
lines with deeper cores than Mm 5394.7 A would show larger 
activity modulation in conflict with the observations. 

Actually, as illustrated in Fig.Q] and demonstrated in 
Figs.HHHl plage thanks its disk-center brightening in any pho- 
tospheric diagnostic not to being hotter at equal height but to 
below- the- surface viewing of hot- wall heat within magnetic con- 
centrations. The SATIRE modeling does not evaluate Wilson de- 
pressions, but as long as the SATIRE models are used in one- 
dimensional radial fashion on their own column mass or op- 
tical depth scale this does not matter. In this sense the Unruh 
plage model does recognize that the local temperature gradient 



around local r = 1 within magnetic concentrations tends to be 
less steep than in the granulation, as in the bottom panels of 
Fig. [71 However, the approximation breaks down for non- vertical 
fluxtube viewing f or which the "Zurich wine-glass" geometry of 
iBiinte et all (1 19931) with slanted rays passing through the glasses 
was a much more realistic description, and it fails for limbward 
faculae because slanted viewing of hot granule innards through 
empty fluxtubes (Fig.[T]) should not be described by the vertical 
temperature stratification within magnetic concentrations. 

Plage and faculae were much better treated by the older 
Solanki-style fluxtube models which diverge with depth instead 
of with height between magnetic and non-magnetic (compare 
the PL A and Unruh plage models in Fig. [2 and PL A with the 
MURaM stratification in Fig. [7]). Obviously, outward divergence 
supplies a zero-order approximation to increasing facular con- 
trast in limbward viewing, but no more than that and inherently 
wrong. 

The same criticisms apply to the similar photospheric feature 
modeling through outward-diverging temperature stratifications 
by e.g., Fontenla et al. (1993, 2006), who add more ad-hoc ad- 
justment parameters in the form of a deep-seated ch romospheric 
tempe rature rise, comparable to the ones invoked bvlDovle et al.l 
(120011) . that was rightfully removed by lUnruhet al.l (1 1999b in 
their plage model - the chromosphere has nothing to do with 
network and plage visibility in photospheric diagnostics. 

Nevertheles s , the success of the SATIRE modeling by 
IDanilovic et al.l (120071) implies that, given any trick to make a 
single magnetic concentration brighter in Mn 1 5394.7 A than in 
Fe 1 5395.2 A by the amount given by SATIRE'S two-model di- 
vergence, that trick will reproduce Livingston's data similarly. 
Thus, although we have not performed any full-disk modeling, 
our trick is likely to reproduce these data too. Our trick entails 
better understanding of why Mm lines activity-brighten more 
than other lines: the latter brighten more in normal granulation. 

5. Conclusion 

The explanation of the activity sensitivity of Mm 5394.7 A con- 
cerns deep-photosphere line formation only. Intergranular mag- 
netic concentrations brighten with respect to field-free intergran- 
ular lanes in any photospheric diagnostic through deep radia- 
tion escape sampling relatively high and flat-gradient temper- 
atures (Figs. [5] and 0). For normal, narrow photospheric lines 
this brightening has less effect in full-disk averaging through 
their loss of line depth in normal granulation (Fig. [7]) that was 
traditionally mimicked by applying micro- and macroturbu- 
lent smearing (Fig. [3]). The flat-bottomed profile which Mm 
5394.7 A possesses thanks to its hyperfine structure (Fig. [3]) 
makes this line much less susceptible to granular Doppler smear- 
ing and thermal broadening so that it weakens less in normal 
granulation (Fig. [7]) and so displays larger mean brightness con- 
trast between quiet and magnetic areas (Fig.|U). 

The Mn 1 5394.7 A line is therefore an unsigned proxy mag- 
netometer sensitive to the magnetic concentrations that consti- 
tute on-disk network and plage and near-limb faculae, through 
Wilson-depression viewing of subsurface bright-wall heat near 
disk center and through slanted facular viewing into hot gran- 
ules near the limb (Fig.Q]). In solar irradiance monitoring Mm 
5394.7 A tracks well with the Ca 11 H & K and Mg 11 h & k core 
brightnesses because these also respond to magnetic concen- 
trations, although through unidentified magnetic chromosphere 
heating that does not affect the Mn 1 lines, neither directly nor 
through interlocking to Mg 11 h & k. 
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As a proxy magnetometer Mm 5394.7 A is similar to the G 
band in which the contrast enhancement arises from the gen- 
eral addition of CH line opacity and local reduction of that 
through dissociation in magnetic concentrations, and to the ex- 
tended blue wings of strong lines in which the contrast enhance- 
ment arises from the general addition of line opacity with reduc- 
tion of that through lesser damping in magneti c concentrations 
plus D oppler flattening of the granular contrast dLeenaarts et alJ 
I2006bl) . These are all sufficiently wide in wavelength to not suf- 
fer from the granular Doppler smearing that spoils the contrast 
for the centers of narrow lines. The G band is the most useful of 
these proxies by being a wide-band spectral feature in the blue. 
Qua contrast, the blue wing of Ha is probably the best of all 
dLeenaarts et al.ll2006al) . but not for full-disk irradiance monitor- 
ing since its photospheric magnetic-concentration brightening is 
sometimes obscured by overlying dark blue-shifted and/or heat- 
widened chromo spheric fibril absorption. 

We conclude that the principal usefulness of photospheric 
Mn i lines lies not in their unusual activity sensitivity but in their 
hyperfine- structured richness as weak-field diagnostic in full- 
Stok es polarimetry with high angu lar resolution and sensitivity 
(e.g. JSanchez Almeida et al.ll2008l) . 
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